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ABSTRACT 

We  are  developing  a  laterally  variable  velocity  model  of  the  crust  and  upper  mantle  across  Eurasia  and  North  Africa 
to  reduce  event  location  error  by  improving  regional  travel  time  prediction  accuracy.  The  model  includes  both  P  and 
S  velocities,  and  we  describe  methods  to  compute  travel  times  for  Pn,  Sn,  Pg,  and  Lg  phases.  For  crustal  phases  Pg 
and  Lg  we  assume  that  the  waves  travel  laterally  at  mid-crustal  depth,  with  added  ray  segments  from  the 
station/event  to  the  mid  crustal  layer.  Our  work  on  Pn  and  Sn  travel  times  extends  the  methods  described  by  Zhao 
and  Xie  (1993).  With  consideration  for  a  continent-scale  model  and  application  to  seismic  location,  we  extend  the 
model  parameterization  of  Zhao  and  Xie  (1993)  by  allowing  the  upper-mantle  velocity  gradient  to  vary  laterally. 

This  extension  is  needed  to  accommodate  the  large  variation  in  gradient  that  is  known  to  exist  across  Eurasia  and 
North  Africa.  Further,  we  extend  the  linear  travel-time  calculation  method  to  mantle-depth  events,  which  is  needed 
for  seismic  locators  that  test  many  epicenters  and  depths.  Using  these  methods,  Sandia  National  Laboratory  has 
developed  a  code  to  compute  regional  travel  times  on-the-fly  from  the  velocity  model  in  milliseconds,  forming  the 
basis  of  a  flexible  travel-time  facility  that  may  be  implemented  in  an  interactive  locator. 

We  use  a  tomographic  technique  to  improve  upon  a  laterally  variable  starting  velocity  model  that  is  based  on 
Lawrence  Livermore  and  Los  Alamos  National  Laboratory  model  compilation  efforts.  Our  tomographic  data  set 
consists  of  approximately  20  million  regional  arrivals  from  events  that  meet  the  ground  truth  (GT)  criteria  of  Bondar 
et  al.  (2004)  and  other  non-seismic  constraints.  Each  datum  is  tested  to  meet  strict  quality  control  standards  that 
include  comparison  with  established  distance-dependent  travel  time  residual  populations  relative  to  the  model 
iasp91.  In  addition  to  bulletin  measurements,  nearly  50  thousand  arrival  measurements  were  made  at  the  national 
laboratories.  The  tomographic  method  adjusts  Pn  velocity,  mantle  gradient,  and  a  node-specific  crustal  slowness 
correction  for  optimized  travel-time  prediction. 
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OBIECTIVES 

This  project  produces  a  laterally  variable  velocity  model  of  the  crust  and  upper  mantle  that  is  specifically  designed 
for  use  in  routine  seismic  location.  At  this  time  the  Seismic  Location  Baseline  Model  (SLBM)  is  focused  on 
travel-time  prediction  at  local  and  regional  distances.  Therefore,  ray  paths  are  wholly  within  the  crust  and  upper 
mantle.  Like  any  travel-time  prediction  method  used  in  a  location  algorithm  the  SLMB  must  return 

1 .  An  accurate  travel-time  prediction  and 

2.  An  uncertainty  estimate  of  the  travel-time  prediction  error 

Because  the  SLBM  is  meant  for  use  in  routine  location  algorithms  where  networks  can  be  dynamic  and  pre- 
computation  of  travel  times  for  all  available  data  may  not  be  possible,  the  SLBM  must  also 

3.  Compute  the  travel  time  on-the-fly  given  regional-  or  local-distance  station/event  coordinates  and 

4.  Return  the  travel  time  in  milliseconds,  thus  enabling  the  estimation  of  a  location  in  a  few  seconds 

Further,  we  aim  to  improve  a  starting  model  that  is  based  on  a  geophysical  compilation.  The  improvement  will  be 
achieved  using  a  ground-truth  data  set  and  a  tomographic  technique  that  is  tailored  to  optimize  model  parameters 
important  to  seismic  location. 

RESEARCH  ACCOMPLISHED 

We  meet  the  objectives  outlined  above  by  adapting  several  approaches  for  model  parameterization  and  travel-time 
calculation  into  one  package  for  computation  of  regional-  and  local-distance  travel  times  for  Pn,  Sn,  Pg,  and  Lg 
phases. 

A  challenge  for  this  project  is  developing  a  model  parameterization  that  enables  fast  and  accurate  prediction  of  each 
local/regional  phase  at  all  applicable  distances.  Well-established  methods  can  be  used  to  compute  regional  Pn  and 
Sn  travel  times  (e.g.,  Hearn,  1984),  but  the  accuracy  of  these  methods  degrades  at  far  regional  (>1000  km)  distance 
(e.g.,  Hearn  et  al.,  2004).  To  more  accurately  predict  Pn  and  Sn  at  far-regional  distances,  Zhao  (1993)  and  Zhao  and 
Xie  (1993)  approximate  upper-mantle  structure  with  a  linear  gradient,  resulting  in  a  simple  expression  for 
calculating  travel  time.  For  Pn  and  Sn,  we  adapt  the  Zhao  and  Xie  (1993)  approach  for  application  to  seismic 
location.  For  Pg  and  Lg  (Sg)  at  local  distances,  we  extract  a  vertical  cross  section  from  the  model  and  use 
2-dimensional  ray  tracing  to  compute  the  travel  time.  At  regional  distance,  we  approximate  the  crustal  waveguide 
(Pg  and  Lg  phases)  with  a  laterally  variable  velocity  layer,  while  accounting  for  propagation  to/from  the 
station/event  using  ray  tracing. 

In  each  instance  we  use  tomographic  methods  to  improve  travel-time  prediction  of  the  model.  Lawrence  Livermore 
and  Los  Alamos  National  Laboratories  (LLNL  and  LANL,  respectively)  are  developing  a  joint  ground  truth  data  set 
for  tomography  that  includes  stringent  quality  control  measures  for  arrival-time  measurements.  We  adapt  Pn 
tomographic  methods  (e.g.,  Hearn  2004;  Zhao  and  Xie,  1993,  Phillips  et  al.,  2007)  to  improve  travel-time 
predictions  of  each  regional  phase. 

M  odel  Parameterization 

We  combine  the  laterally  variable  layer  approach  of  Pasyanos  et  al.  (2004)  with  the  linear  mantle  gradient  of  Zhao 
and  Xie  (1993).  Layer  definitions  are  specified  in  Table  1.  Note  that  the  thickness  of  some  layers  may  be  zero.  For 
instance,  on  the  continents  the  depth  of  the  water  layer  coincides  with  the  depth  of  the  model  layer  exposed  at  the 
surface.  Velocity  vs.  depth  profiles  are  defined  at  nodes,  and  the  profiles  at  the  nodes  are  interpolated  using  an 
efficient  code  developed  at  Sandia  National  Laboratory  (SNL)  to  determine  velocity  at  any  arbitrary  location 
(lat,lon,depth).  SNL  has  also  developed  a  tessellation  node  structure  on  a  spheroid  with  node  spacing  of 
approximately  1°  (Figure  1).  At  present,  the  model  development  domain  is  Eurasia  and  North  Africa,  and  nodes 
inside  that  domain  capture  the  effects  of  3 -dimensional  structure  on  travel  times.  Outside  of  the  development  domain 
nodes  are  set  to  a  default  velocity  profile  based  on  iasp91  (Kennett  and  Engdahl,  1991).  This  parameterization 
provides  a  seamless  and  extensible  model.  Expansion  beyond  Eurasia  and  North  Africa  does  not  require  a  change  in 
the  model  parameterization  itself,  only  modification  of  the  velocity  structure  at  previously  defined  nodes.  Further, 
SNL  has  incorporated  the  GRS80  ellipsoid  in  to  the  model,  eliminating  the  need  for  the  conventional  ellipticity 
correction  to  travel-time  predictions. 
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Figure  1.  SL  BM  global  parameterization,  a)  An  example  tessellation  with  approximately  1°  grid  spacing. 

Color  is  based  on  approximate  M  oho  depth,  b)  An  example  velocity/depth  profile  as  defined  at  each 
node.  T  he  mantle  portion  of  the  profile  is  specified  by  the  velocity  at  the  crust/mantle  interface  and 
a  linear  gradient. 

T  able  1.  M  odel  entities  and  associated  parameters.  Depths  and  velocities  at  each  node  are  interpolated  to 

define  a  3-dimensional  crustal  model  that  overlies  a  laterally  variable  model  in  the  shallow  mantle. 
Notethat  the  mid-crustal  layer  isdistinct,  in  that  the  velocities  are  specified  individually  for  each 
phase.  For  Pn  and  Sn  the  mid-crustal  layer  is  used  in  conjunction  with  all  other  crustal  layers  to 
compute  travel  times  for  rays  that  travel  steeply  through  the  crust.  For  Pgand  Lg  distinct  middle 
crust  velocities  are  used  to  model  the  horizontal  slowness  of  the  regional  phases  that  are  trapped  in 
the  crust. 


Model  entity 

Layer 

Depth 

P-wave  velocity 

S-wave  velocity 

P-wave 

gradient 

S-wave 

gradient 

Water 

Yes 

Yes 

Yes 

Sediment  1 

Yes 

Yes 

Yes 

Sediment  2 

Yes 

Yes 

Yes 

Sediment  3 

Yes 

Yes 

Yes 

Upper  Crust 

Yes 

Yes 

Yes 

Middle  Crust 

Yes 

Independent  for  Pn,  Pg 

Independent  for  Sn,  Lg 

Lower  Crust 

Yes 

Yes 

Yes 

Moho 

Yes 

Yes 

Yes 

Mantle  Gradient 

Yes 

Yes 

Travel-Time  Calculation,  Pn  and  Sn 

The  travel-time  calculation  is  based  on  the  method  described  in  Zhao  (1993)  and  Zhao  and  Xie  (1993).  This 
calculation  is  similar  to  the  widely  used  approach  of  Hearn  (1984),  with  an  additional  term  (y)  introduced  to  account 
for  diving  rays  that  may  occur  due  to  a  positive  velocity  gradient  with  depth  and  Earth  sphericity. 
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The  travel-time  calculation  is 

TT  =  Y.diSi+ CJ(+ 13+ Y  [1] 

i=\ 

where  d  and  ^  are  the  distance  and  slowness  (taken  as  1/Moho  Velocity)  in  each  of  the  i  segments  comprising  the 
great-circle  path  between  Moho  pierce  points  near  the  station/event,  a  and  (3  are  the  crustal  travel  times  at  the  source 
and  receiver,  and  y  is  a  term  that  accounts  for  the  effect  of  both  mantle  velocity  gradient  and  earth  sphericity. 

We  define  a  as 


M 
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where  v  and  r  are  the  velocity  and  layer  radius  of  the  M  crustal  layers  from  the  event  to  the  Moho,  and  p  is  the  ray 
parameter  (p=l/v,  v  evaluated  at  the  ray  bottoming  depth). 


We  similarly  define  (3  as 


N 

k=\ 


7  2  P  1  /  2  P 

/Vk  \  /Vk 


[3] 


where  v  and  r  are  defined  as  above  for  the  L  crustal  layers  from  the  station  to  the  Moho.  The  same  p  is  used  in  both 
Eqns  [2]  and  [3].  Because  the  ray  bottoming  depth  is  a  function  of  the  pierce  point,/?  is  determined  through  an 
efficient,  iterative  process. 


Per  Zhao  and  Xie  (1993), 


Y=  — 


c^Xm^ 


[4] 


24  Fo 

where  is  the  horizontal  distance  traveled  in  the  mantle,  c  is  a  velocity  gradient  in  the  mantle  that  is  normalized  by 
the  velocity  at  the  crust  mantle  boundary  plus  an  additional  term  to  account  for  Earth  sphericity  (Helmberger,  1973), 
and  Vo  is  a  regional  average  of  Pn  velocity  over  the  entire  study  area. 


We  introduce  spatially  varying  c  into  the  model  (Phillips  et  al.,  2007),  and  we  calculate  y  by  averaging  c  along  each 
ray.  Vq  remains  an  average  Pn  velocity  over  the  whole  model,  which  allows  us  to  take  advantage  of  linear 
tomographic  inversion  methods  (see  below).  Tests  suggest  that  the  approximation  to  Vq  introduces  negligible 
travel-time  error  given  Pn  velocities  ranging  from  7.5  km/s  to  8.3  km/s. 


The  Zhao  (1993)  method  is  applicable  to  events  in  the  crust,  making  the  approach  well  suited  to  nuclear  explosion 
monitoring.  However,  seismic  location  algorithms  may  explore  the  possibility  that  an  event  occurred  in  the  mantle, 
necessitating  a  consistent  method  of  travel-time  predictions  for  mantle  events.  The  following  extends  the  travel-time 
method  to  events  in  the  shallow  mantle,  with  the  condition  that  c^h^«l  (h  is  the  bottoming  depth  of  the  ray) 


TT  —  (X+  tyr 


[5] 


where  a  is  the  crustal  travel  time  from  the  Moho  to  the  station  (as  defined  in  [2]),  and  4  is  the  travel  time  in  the 
mantle.  Figure  2  shows  the  geometry  and  defines  many  of  the  variables  used  in  the  following  equation. 
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[6] 


If  the  ray  leaves  the  event  upwards,  then  the  second  term  is  subtracted.  If  the  ray  leaves  the  event  downwards,  then 
the  second  term  is  added,  tuoko  is  the  travel  time  for  a  ray  traversing  the  Moho  from  the  event  to  the  point  where  the 
ray  enters  the  crust  and  propagates  to  the  station,  is  the  horizontal  distance  as  measured  at  Moho  radius  by  a  ray 
that  starts  at  the  Moho  then  travels  downward  passing  through  the  event  and  continuing  to  the  station,  x^  is  similar  to 
x^,  but  the  horizontal  distance  is  measured  at  the  radius  of  the  event,  d  is  the  horizontal  distance  traveled  in  the 
mantle  from  the  event  to  the  Moho  pierce  point  below  the  station,  as  measured  at  Moho  radius,  is  the  mantle 
velocity  gradient  normalized  by  average  Moho  velocity,  with  the  addition  of  a  term  to  account  for  earth  sphericity 
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(Helmberger,  1973).  is  the  depth  of  the  event  below  the  Moho.  Voz  is  the  average  model  velocity  at  the  depth  of 
the  event. 


Upgoing  CaseM 


Downgoing  Case 


Event 


F  igure  2.  G  eometry  and  variable  definition  extending  the  Z  hao  and  X  ie  (1993)  formulation  to  events  in  the 
shallow  mantle.  The  red  star  is  the  event  location  and  the  triangle  is  the  station  location. 

Travel-Time  Calculation,  Pg  and  Lg 

Both  Pg  and  Lg  phases  are  trapped  in  the  crust,  and  both  phases  exhibit  complex  waveforms  that  require  hundreds  or 
thousands  of  rays  to  model.  Further,  the  first  arriving  ray  with  sufficient  energy  to  be  observed  is  dependent  on 
geologic  structure  and  station/event  distance.  Therefore,  it  is  difficult  to  physically  model  the  observed  travel  time 
using  conventional  ray  techniques.  Empirically,  Pg  and  Lg  travel  at  horizontal  velocities  of  approximately  6.0  km/s 
and  3.2  km/s  (respectively),  which  is  suggestive  of  propagation  in  the  middle  crust.  Further,  event  depth  can  impart 
a  static  travel-time  delay,  suggesting  a  component  of  propagation  from  the  event  to  the  middle  crust  then  up  to  the 
station.  We  capture  this  travel-time  behavior  with  a  simple  approximation  whereby: 

TT^YdiSi  +  a+P  [7] 

where  di  is  the  distance  traveled  in  the  middle  crust  in  each  of  N  ray  segments  and  a  and  (3  propagate  the  phase  to 
and  from  the  middle  crust,  respectively.  When  the  source  is  above  the  middle  crust,  then  the  calculation  is  almost  the 
same  as  [1],  but  the  correction  for  a  diving  ray  is  not  used.  When  the  source  is  below  the  mid-crustal  layer,  we 
assume  that  the  ray  travels  horizontally  until  Earth  sphericity  causes  the  ray  to  intersect  the  mid-crustal  layer.  While 
this  approach  by  no  means  captures  the  physical  complexity  of  Pg  and  Lg  wave  propagation,  we  find  that  it  is 
suitable  for  estimating  travel  time. 
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Computational  Efficiency 

Travel-time  facilities  used  in  routine  location  algorithms  must  be  computationally  efficient.  Analysts  often  iterate  on 
arrival-time  measurements,  first  locating  the  event  with  clear  arrivals,  then  adding  additional  arrivals  (or  adjusting 
previous  measurements)  based  on  a  preliminary  location.  The  demands  of  routine  seismic  analysis  preclude  waiting 
more  than  a  second  or  two  for  the  calculation  of  a  location  before  continuing  work.  Considering  that  a  location 
requires  the  computation  of  travel  times  for  many  phases,  over  many  iterations,  often  with  numerical  calculation  of 
derivatives  (at  least  4  travel-time  calculations  per  observation  per  iteration),  the  time  to  compute  an  individual  travel 
time  must  be  negligible  (milliseconds). 

The  algebraic  form  of  the  travel-time  formulas  that  are  specified  above  consume  negligible  computer  time.  The 
primary  challenge  is  fast  extraction  of  the  cross  sectional  profile  needed  for  the  travel-time  calculation  from  the 
laterally  variable  model.  The  tessellation  model  parameterization  developed  at  SNL  takes  advantage  of  mature 
algorithms  to  determine  which  triangle  any  given  point  lies  within,  which  nodes  comprise  the  triangle,  and  node 
weights  used  for  interpolation.  Figure  3  shows  an  example  of  a  cross  section  and  the  Pn  ray  that  is  used  to  compute 
the  travel  time. 

Test  results  show  a  linear  increase  in  computer  time  as  a  function  of  station/event  distance.  In  the  distance  range 
from  -200  km  to  -2200  km,  Pn/Sn  travel  times  require  -1  milliseconds  to  2  milliseconds  to  compute  on  a  desktop 
computer  with  -1.5Ghz  clock  speed.  Pg/Lg  times  are  faster. 
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Figures.  Cross  section  extracted  from  theiateraiiy  variabieSLBM  modei.  The  components  of  the  Pn/Sn 

travei-time  calculation  are  also  shown.  The  blue,  red,  green,  and  cyan  colors  correspond  to  the  first, 
second,  third,  and  fourth  terms  of  Eqn  1. 

T  ravel-Time  Accuracy 

We  make  use  of  many  approximations  to  meet  the  challenge  of  fast,  on-the-fly  travel-time  calculation.  We  test  the 
accuracy  of  these  calculations  by  comparing  travel  times  computed  using  the  methods  described  above  (SLBM)  with 
travel  times  computed  using  a  fully  3 -dimensional  calculation  described  in  Flanagan  et  al.  (2007).  The  Flanagan  et 
al.  (2007)  study  uses  the  WENAl.O,  3-dimensional  model  of  Pasyanos  et  al.  (2004).  For  testing  purposes,  we  use  the 
WENAl.O  crustal  model  and  mantle  velocity  at  the  Moho  discontinuity.  We  then  compute  the  linear  mantle  gradient 
using  the  velocity  at  the  Moho  and  the  velocity  at  130  km  depth. 

Figure  4  shows  representative  results  of  our  tests.  Figure  4a,b  shows  the  difference  in  Pn  travel  times  computed 
using  the  SLBM  method  and  using  3 -dimensional  finite  difference,  respectively.  In  both  cases  the  comparisons  are 
relative  to  the  iasp91  model,  and  Pn  travel  times  are  computed  using  the  TauP  toolkit  of  Crotwell  et  al.  (1999).  The 
general  features  of  Figure  4a  and  4b  are  similar,  with  deviations  from  iasp91  travel  times  ranging  from 
approximately  2  seconds  late  to  6  seconds  early.  Figure  4c  shows  the  difference  between  the  SLBM  and  finite- 
difference  calculations,  with  deviations  between  the  two  methods  generally  between  plus  and  minus  1  second.  Errors 
in  the  finite  difference  calculations  themselves  are  reported  to  be  approximately  plus  and  minus  0.5  seconds 
(Flanagan  et  al.,  2007).  These  tests  suggest  that  the  efficient  SLBM  computational  errors  are  slightly  larger  than  the 
computationally  demanding  finite-differences  calculations,  and  SLBM  errors  are  relatively  small  compared  to 
iasp91  travel-time  anomalies. 
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Figure4.  Comparision  of  Pn  travel  times  using  SLBM  and  3D  finite  difference  calculations  for  station  BGCA 
in  central  Africa.  For  SLBM  and  3D  finite-difference  calculations  the  W  ENA  1.0  model  of  Pasyanos 
etal.  (2004)  is  used.  Circles  around  the  station  are  5°  increments.  See  text  details. 

T  omography 

Test  results  presented  above  suggest  that  the  SLBM  travel-time  method  is  in  good  agreement  with  3 -dimensional 
calculations.  The  WENAl.O  model  used  in  the  tests  above  is  shown  to  improve  travel-time  prediction  and  location 
accuracy  relative  to  the  iasp91  default  model  (Flanagan  et  ah,  2007).  We  will  further  improve  SLBM  travel-time 
accuracy  using  tomographic  methods. 

A  data  set  with  small  errors  in  event  location  and  arrival-time  measurements  is  critical  to  tomographic  studies.  Data 
coverage  is  also  critically  important.  LLNL  and  LANL  have  combined  ground-truth  data  sets  for  this  study.  Both 
national  laboratories  contribute  global,  regional,  and  local  bulletins  (some  not  widely  available),  as  well  as  tens  of 
thousands  of  arrival-time  measurements  made  at  the  national  laboratories.  All  event  locations  are  evaluated  against 
Bondar  et  al.  (2004)  epicenter  accuracy  criteria,  and  all  picks  are  evaluated  against  an  error  budget  that  accounts  for 
event  mislocation,  iasp9 1  prediction  error,  and  arrival-time  measurement  error.  Observations  outside  of  the  99% 
confidence  bounds  for  total  error  are  removed.  Figure  5  shows  the  extensive  data  coverage  for  Pn,  Sn,  Pg,  and  Lg 
throughout  Furasia. 

We  are  using  the  SLBM  travel-time  calculator  in  tomographic  inversion  programs.  The  general  form  of  the 
tomographic  inversion  for  Pn/Sn  is  as  follows: 
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A  significant  difference  between  the  formulation  presented  here  and  more  typical  Pn-tomography  formulations  is  the 
introduction  of  a  scalar  value  that  adjusts  the  slowness  of  the  crustal  stack,  as  opposed  to  a  static  time-term  to 
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account  for  crustal  travel-time  delays.  Adjusting  the  crustal  slowness  produces  a  model  that  is  better  suited  to 
account  for  the  effects  of  event  depth  on  predicted  travel  times  and  therefore  better  for  estimating  event  depth. 

The  tomographic  free  parameters  are  more  easily  identified  in  matrix  form: 
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Regularization 


where 


T  =  travel  time, 

s  =  Pn  slowness, 

X  =  Pn  distance  (or  weight), 

c  =  normalized  velocity  gradient  v=vo(l  +cz), 

Xm  =  length  of  Pn  path, 

Vo  =  average  Pn  velocity, 

V  =  velocity  of  a  crustal  layer, 

k  =  index  on  K paths  (travel-time  observations), 

i  =  index  on  N  model  nodes  (mantle  path), 

j  =  index  on  N  model  nodes  (crustal  path), 

p  =  index  on  Q  crustal  layers,  and 

a  =  scalar  adjustment  to  the  crustal 

slowness  stack  at  each  node. 


CONCLUSIONSAND  RECOMMENDATIONS 


We  describe  the  progress  of  the  SLBM  project  to  date.  This  project  is  distinct  because  it  tailors  the  travel-time 
prediction  algorithm  and  tomography  results  for  use  in  routine  seismic  location  algorithms.  Emphasis  is  placed  on 
travel-time  prediction  accuracy  and  computational  efficiency  of  regional  phases. 

The  tessellation  model  parameterization  provides  seamless  global  coverage.  The  use  of  a  tessellation  approach  also 
allows  fast  interpolation  of  model  parameters  to  extract  the  great-circle  cross  section  of  velocity  structure  that  is 
needed  to  compute  regional  travel  times.  The  current  focus  of  the  SLBM  effort  is  Eurasia  and  North  Africa,  and 
model  nodes  outside  of  that  area  are  set  to  a  default  velocity  structure  based  on  the  iasp91  model 
(Kennett  and  Engdahl,  1991).  We  note  that  the  model  is  extensible,  and  a  global  calibration  effort  would  entail 
updating  node-centered  velocity  profiles  (by  whatever  means),  which  does  not  require  updates  to  the  model 
tessellation  or  travel-time  codes. 
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F  igure  5.  T omographic  data  set.  Purple  triangles  are  stations  with  at  least  1  regional-distance  arrival-time 
measurement.  G  round-truth  location  and  accuracy  (as  defined  in  Bondar  etal.,  2004)  are  shown  in 
the  legend  in  the  lower  left.  Regional  ray  coverage  (shown  by  gray  lines)  is  excellent  throughout  the 
study  area,  with  the  exception  of  North  Africa. 
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We  make  use  of  several  approximations  that  result  in  a  relatively  simple  algebraic  form  for  travel-time  calculations. 
We  have  tested  the  accuracy  of  the  approximate  methods  against  full,  3 -dimensional  finite-difference  calculations. 
The  differences  between  the  approximate  methods  and  full,  3 -dimensional  finite-difference  methods  are  estimated  to 
be  less  than  10%  of  errors  observed  by  using  a  simple  iasp91  background  model.  Therefore,  we  conclude  that  the 
approach  developed  here  is  a  significant  improvement  in  routine  location  practice. 

We  have  developed  an  extensive,  quality-controlled,  ground-truth  data  set  across  Eurasia  and  North  Africa.  This 
data  set  will  be  used  in  tomographic  codes  that  employ  the  SLBM  travel-time  calculator  to  further  optimize  the 
model  parameters  that  are  most  important  to  seismic  location. 
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